function y = f(x)
    y = 1/sqrt(1 - (x ^2 * sin(x)^2));
endfunction
